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ABSTRACT 

We perform a numerical simulation of the corona of the young, rapidly rotating KO dwarf AB 
Doradus using a global MHD model. The model is driven by a surface map of the radial magnetic 
field constructed using Zeeman-Doppler Imaging. We find that the global structure of the stellar 
corona is dominated by strong azimuthal tangling of the magnetic field due to the rapid rotation. The 
MHD solution enables us to calculate realistic Alfven surfaces and we can therefore estimate the stellar 
mass loss rate and angular momentum loss rate without making undue theoretical simplifications. We 
consider three cases, parameterized by the base density of the corona, that span the range of possible 
solutions for the system. We find that overall, the mass and angular-momentum loss rates are higher 
than in the solar case; the mass loss rates are 10 to 500 times higher, and the angular momentum 
loss rate can be up to 3 x lO"* higher than present day solar values. Our simulations show that this 
model can be use to constrain the wide parameter space of stellar systems. It also shows that an MHD 
approach can provide more information about the physical system over the commonly used potential 
field extrapolation. 
Subject headings: stars: coronae - stars: activity - stars: magnetic field 



1. INTRODUCTION 

Rapidly rotating young stars are important for the 
study of stellar activity in two aspects. First, these stars 
provide an "enhanced" picture of the fundamental stel- 
lar magnetic activity also seen in less active stars (such 
as our Sun) and therefore help to better understand the 
causes and consequences of stellar activity. Second, these 
rapidly rotating systems provide an opportunity to ob- 
serve the repeating rotationally modulated signatures of 
activity on relatively short time scales. 

An example of a well observed young, active star is 
AB Doradus (HD 36705, AB Do r hereafter), a KO dwarf 
with an age of about 75 Myr fZuckermanet al. 20041 
Luhman et al. 2005; Nielsen ct al. 2005; Lopez-Santiago 
et al.|2006| |Janson et al.||2007|) . AB Dor spi ns with a ro 



tation period of J^ = 0.5 dajys (|Pakull||1981 ), has a mass 
of about M* = O.76M0 CCuirado et al."1997"), and ra- 
dius of about R^, = O.86i?0 (Ma ggio ct al. 2000 ) (slightly 
smaller than the solar radius, Rq, and solar mass Mq). 
AB Dor has been well o bserved over mos t of the e lectro- 
magnetic spectrum (e.g. Lim et al.||1994 Grothuc s et al. 
TTuT^ io„i,.„:xx ^^ -1 " 19971 iviihu et ai."1998; Cu tispoto 

and its 



1997 



1998 
activi 



Schmitt et al 



Jarvinen et al. 



]^IViiliu et al.' 
5]n3udding ct al. 



;y has been stuc 
Maggio etal.|2000; Giid el et al..2001^ 



2005J pudding ct al. 2009) 

led m detail in the X-ray band (e.g. 

Sanz-Forcada et al.l 



2003; Garcia- Alvarez e t al.||20051 JHussain et al.| |20051 



Matranga et al.||2005| |Hussain et al. |2007p . Most impor- 
tantly tor the work in nand, AB Uor nas been extensively 



tantly tor the work in nand, AB Uor 

observed using Zeeman-Doppler I maging (ZDI) ( [Donati 

& Semel||1990 ponati et al.||2009|). This method enables 



some inference ot the vector magnetic field on the stellar 
surface based on the polarization of emitted light by the 
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field and the modulation in Zeeman splitting as the star 
rotates. For rapidly rotating systems such as AB Dor, it 
is possible to use this method to reconstruct a large-scale 
map of the stellar surface magnetic topology, analogous 
to the solar synoptic "magnetograms" (27-day~average 
surface field distribution), albeit with lower spatial reso 
lution. The ZDI technique has been used to m ap the sur 



face magnetic field of AB Dor ( Donati & Collier Cameron 
1997| [Donati et al.||1999[ IHussaiiTera l. 2002), as well to 



study its stellar cyclic activity and surface differential ro 
tation (Donati & Collier Cameron 1997; Collier Cameron| 
fc Donati||2002j |Pointer ct al.||2002, |Jefters et al.|.2 00 
I'he stellar cycle ot AB Dor has been also studie J 



onl 



usin g 
Doppler Imaging observations by [Jarvinen et alT (2005). 
It is commonly assumed that, like the solar corona, 
stellar coronae are dominated by their magnetic fields, 
so that the magnetic pressure, Pb = B'^ /Sir, is much 
greater than the thermal pressure, Pth = nkT, and that 
the plasma /3 parameter, f3 = Pth/Ps, is much smaller 
than 1. In this case, the magnetic field can be assumed 
to be potential (i.e. there are no forces or currents acting 
on it) and it can be described as a gradient of a scalar po- 
tential. Under these assumptions, the three-dimensional 
distribution of the magnetic field can be obtained by solv- 
ing Laplace's equation for the scalar potential, where the 
surface field maps are used as the inner boundary condi- 
tion and the outer boundary condition assumes a purely 
radial field at a certain height above the surface (known 
as the "source surface"). This technique of extrapolat- 
ing the coronal magnetic field is known as the "potential 



field" method (Altschuler & Newkirk 1969 Altschuler 



et al.|197"7 1 . It has been used extensively m solar studies 



and more recently to extrapola te the coronal magnetic 
fields of stars such as AB Dor ( Donati et al.[[1999[ [Jar 



s 



dine et al.|1999l jHussain et al.||2002jJJardine et aL|2002 



Mclvor et al.[[2003[ ). In particular, [Hussainet"ar[|2007l 
have reconstructed the X-ray corona of AB Dor based on 
the potential field extrapolation, combined with different 



X-ray models for coronal loops. 

The potential field extrapolation is a useful tool to ob- 
tain a first order approximation of the large-scale struc- 
ture of a stellar corona bas ed on its surface ma gnetic 
map. The approach taken by'Hussain et al. (2007) is jus- 
tified due to the fact that the major part ot the coronal 
X-ray emission is expected to originate from the smaller 
closed loops near the surface, which are usually in a near- 
potential state (in the static case where footpoint mo- 
tions and other short-term motions are not taken into 
account). There are, however, good reasons to attempt 
a more physical approach for describing stellar coronae. 
First, the potential field approximation (by itself) pro- 
vides information only about the magnetic field of the 
system, and does not address energy dissipation through 
driving a wind. Second, the location of the source surface 
is not well defined, and third, when considering a com- 
plete description of the physics involved, including con- 
servation of mass, momentum, and energy, one needs to 
take into account the effects of coronal heating and stel- 
lar wind acceleration, the stretching of the field lines by 
the highly conductive coronal plasma to a non-potential 
state, as well as the effects of rapid rotation in stars like 
AB Dor. 



Here we extend the work of Hussain et al. ( 2007 ) and 
present a complete three-dimensional MagnetoHydroDy- 
namic (MHD) simulation of the corona of AB Dor based 
on its observed surface magnetic field distribution. For 
our simulation, we use a global MHD model developed for 
the solar corona, which provides a self-consistent stellar 
wind solution driven by surface magnetic field maps. The 
end result is a steady state, MHD, non-potential solution 
of the corona and wind of AB Dor, which includes the 
distribution of the complete set of physical parameters 
in the simulation domain. This more complete solution 
provides a better understanding of the large-scale coro- 
nal structure. We highlight the differences between the 
coronae of young stars like AB Dor and the solar corona 
due to rapid rotation of the former. We also provide re- 
alistic calculations of the possible mass loss rates for AB 
Dor, parameterized by the coronal base density. 

We present the numerical model and the observational 
constraints used in the simulation in Section [2] The re- 
sults are presented in Section [3] and the main findings 
are discussed in Section |4l We conclude this work in 
Section [5] 

2. NUMERICAL SIMULATION 

The simulatio n of AB Dor is done u sing the solar 
corona model by [Cohen et al. ( 2007 2008 1 , which is part 
of the Space Weather Modeling Framework (SWMF) 
(|Toth et al.|[2005) and is ba sed on the g eneric MHD 
BATS-R-US model by ( jPoweU et al.||1999[ ). The model 
is driven by surface magnetic field maps, and the initial 
condition for the magnetic field, as well as the volumetric 
energy input for the stellar wind acceleration, is based on 
the distribution of the potential field. In addition, the 
boundary condition for the surface plasma density, po, 
is scaled with the magnetic field so that the plasma at 
closed field regions is more dense t han in open field re - 
gions, as observed for the solar case ( Phillips et al.|1995 ). 

In the solar case, the source surface is usually set to be 
at r = 2.5i?0. In the case of AB Dor however, the surface 
distribution of the magnetic field contains large regions 



with strong field. Therefore, we expect loops on AB Dor 
to be much larger than solar loops so we choose to set 
the source surface at r = lOi?*. This should not have 
any effect on the non-potential, MHD solution since the 
potential field only serves at the initial condition. The 
MHD solution is mostly affected by the distribution of 
energy deposited into the stellar wind, and this energiza- 
tion is not sensitive to the location of the source surface 
as long as it is set above the the height of the largest 
closed loops. However, setting the source surface below 
the actual size of the loops (at r = 2.5i?j, for example), 
forces more field lines to be open and as a result, each 
plasma cell in the stellar wind is over-energized, resulting 
in solutions with unrealistically fast stellar winds. 

A self-consistent wind acceleration in the code is ob- 
tained by assuming an empirical relation between the 
magnetic flux tube expansion and the terminal stellar 
wind o rigi nating from that flux tube. |Wang fc Sheeley] 



( |1990| and jArge fc Pizzo 
cal formula that relates t 



( 2000 1 have derived an empiri 
le final solar wind distribution, 
Uswi to the flux tube expansion factor, fa- The factor /s 
is the ratio of the magnetic flux of a particular flux tube 
at r = Rgs and at r = i?^,, where Rgg is the height of the 
source surface. The empirical method described above 
predicts the spherical distribution of the solar wind speed 
at r — >■ cx). It is reasonable to assumed that far from the 
Sun (or star), the total energy equals to the bulk kinetic 
energy of the plasma, while on the solar surface, t he 
total energy equals to the enthalpy of the fluid, minus 
the gravitational potential energy (the kinetic energy is 
zero) . By adopting the conservation of total energy along 
a streamline (Bernoulli Integral), we can relate the final 
solar wind speed, Usw and the surface value of polytropic 
index, 70, assuming the boundary conditions for the sur- 
face temperature, Tq, are known: 



2 70 - 1 rup Ri, 



70 



e-U 



with 



hTo 



sw 



GAL 
R. 



(1) 
(2) 

(3) 



Here k^ being the Boltzmann constant, nip the proton 
mass, and G the gravitational constant. 

Close to the Sun, the value of 7 is observed to be about 
unity, (the plasma is hig hly turbulent), whi l e at 1 AU 7 
has a value closer to 1.5 (iTotten et al. 119951 119961). This 



observed modulation in 7 can be related to the power- 
ing of the solar wind, in the manner that the larger the 
gradient in 7 along a flux tube, the faster the wind flows 
along that tube. Based on this assumption, and on the 
relation presented in Eq. [T] it is possible to construct a 
volumetric heating function, _E^(7o, r), in a way that the 
observed volumetric acceleration of the solar wind can 
be recovered. The additional term S^ — j' as 7 — ?► 3/2. 
The model described here constructs the particular 
spatial distribution of E^ based on the input surface 
magnetic map and its potential field. It then solves the 
set of conservation laws for mass, momentum, magnetic 
induction, and energy (the ideal MHD equations): 
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with r = 3/2 until a steady-state stellar wind solution 
is obtained. The stellar input parameters required for 
the model are the boundary value for the density, po, as 
well as the stellar radius, i?*, mass, M^,, and rotation fre- 
quency, ri*. A summary of the stellar parameters of AB 
Dor adopted for the simulation, based on the references 
cited in Section [TJ is provided in Table [T] 

In the simulations presented here, we assume that the 
relation between the flux tube expansion and the ter- 
minal speed obtained from that flux tube is a universal 
process that occurs on AB Dor in a similar manner to 
the Sun. Observations of the corona of AB Dor indicate 
dominant plas ma temperatures peaking in the range 3- 
30 MK (e.g. iSanz-Forcada et al.l 12003 iGarcia-Alvarez 



g. iSanz- 



et al.||2005[ ). I'o is, in principle, the average tempera- 
ture oi the stellar corona. However, we stress that in our 
model. To is essentially a free parameter for the bound- 
ary condition that controls the energization of the stellar 
wind. Further description of the adaptation of a solar 
corona mod el to stellar coronae can be found in (Cohen 
erar][2010| 



Figure"^ shows the the input surface magnetic field 
map adopted for the simulations. This is based on spec- 
tropolarimetric observations obtained in December 2007 
(Hussain et al. submitted to MNRAS), and analyzed 
in a similar manner to the maps described in [Hussain] 
et al.| ( |2002[ |2005[|2007| . The reader is referred to those 
works for further details. Since AB Dor has an incli- 



nation of 60 deg ( Kuerster et al. 1994 ) , the part of the 
stellar surface near the far pole is always hidden from 
view. Consequently, surface maps for AB Dor are in- 
trinsically incomplete. The initial, incomplete map is 
shown in the top left panel of Figure [l] To construct a 
complete map, we enforced hemispherical reflection sym- 
metry on the magnetic field across the equatorial plane. 
Those parts of the southern hemisphere with magnetic 
field magnitude of less than 50 G were assigned magnetic 
field values from the same longitude at the corresponding 
northern hemisphere latitude, but with the opposite po- 
larity. The complete map used in the simulation is shown 
as a longitude-latitude contour map (top-right panel) and 
as spherical plots of the two longitudinal hemispheres, 
colored with contours of the surface magnetic field (lower 
panels). Figure [2] shows the three-dimensional distribu- 
tion of the potential field calculated based on this input 
surface map. It can be seen that the loops extend up to 
the height of the source surface (located at r = lOi?*) 
and that they have no toroidal component. The field is 
fully radial above the source surface as required by the 
analytical solution. 

We caution that there are a number of sources of sys- 
tematic uncertainty present in the simulations. First, 
the ZDI maps are missing part of the stellar surface, and 



we have extrapolated the field, assuming antisymmetry, 
to those non-visible areas. Second, the magnetic field 
maps are of limited resolution (latitudinal resolution of 
3 degrees) and do not resolve the details of the active 
regions. Third, large areas in the maps that appear to 
have strong fields may in fact be dominated by local- 
ized active regions that are smeared out due to the lower 
resolution; such a scenario would lead to a significantly 
different MHD solution. While the ZDI maps do likely 
miss low level structure, the maps should recover the 
strongest field regions, which are likely to dominate over 
the large-scale global models such as those considered 
here. The weaker complex fields would only be interest- 
ing much closer to the surface. 
Observations of AB Dor reveal that the coronal base 



density, np, ranges b etween 10 



10 



10 



12 



( Sanz- 



Forcada et al. 20031. These measurements, however, 
were made based on measurements of strong emission 
lines that are associated with the denser, closed loops. 
Therefore, the density in the "quiet star" (analogous to 
the "quiet Sun") where the wind originates should be 
lower. In order to partially cover this density range, 
we simulate three test cases with different values for the 
coronal base density. "Case A", with rig = 2 • 10^ cm~^, 
which is the value used for simulations of the solar wind, 
"Case B" with no = 10^ cto-^, and "Case C" with 
no = 10""^° cm~^. We simulate the wind and corona us- 
ing a Cartesian box of 30i?*x30i?*x30i?*, in the frame of 
reference rotating with the star (to expedite convergence 
using the local time step algorithm ( ^Cohen et al.^2008) ). 
We use a non-uniform grid with a maximum resolution 
of 2 • 10~^i?Vf prescribed near the surface. The grid is dy- 
namically refined during the simulation so that high res- 
olution is applied at location of magnetic field inversion 
(current sheets). We performed the simulation using the 
PLEIADES super computer at the NASA AMES center. 

3. RESULTS 

The steady state MHD solutions for the three test cases 
are shown in Figure [31 The most notable feature of the 
solutions for all cases is the tangling of the field in the 
azimuthal direction due to the rapid rotation of the star. 
This feature clearly does not appear in the potential field 
solution, for which stellar rotation is not a relevant pa- 
rameter, nor in sim ilar MHD solutions for the Sun (e.g. 
Cohen etal?|[2008| . 

As might be expected, solutions for the three cases are 
qualitatively quite similar, though closer inspection does 
reveal significant differences. It can be seen from the 
middle panel of Figure |3] that the radial wind speed de- 
creases with an increase of the base density. In addition, 
the corona is denser in the solution for Case C as com- 
pared to Case A. While many of the tangled field lines in 
Cases A and B are open due to the strong radial stretch- 
ing by the stellar wind, in Case C, most of the tangled 
field lines are closed. The closed loops in the low corona 
are radially stretched in Cases A and B, while the same 
loops are more potential and less stretched in Case C as 
seen in the bottom panel of Figure [3] 

The interplay between the radial speed and the coronal 
density structure determines the stellar mass loss rate, as 
well as the stellar angular momentu m loss rate to the stel - 
lar wind. We follow the method by Cohen et al. ( [2009 ), 
and calculate these loss rates from the MHD solution. 



This method ex pands the ideahzed approach by Weber| 
& Davis (19671 and uses the fact that the MHD solu- 
tion provides a reahstic, non-ideahzed Alfven surface, at 
which the Alfvenic Mach number, Ma = u/va = 1, 
where va = B / ^/Anp is the Alfven speed. Once the 
Alfven surface has been determined, the loss rates can 
be calculated as: 



M 



pu • daA, 



J 



il^ sin 9 Ta pu- dsLA, 



(5) 



(6) 



where ta is the local radius of the Alfven surface, da^ 
is a surface element, and the integration is done over 
the realistic Alfven surface. It is worth to mentioning 
that the realistic Alfven surface, at which the magnetic 
breaking of the stellar wind takes place, is the actual 
source surface, and it does not have a spherical shape as 
is assumed to have in the common use of the potential 
field approximation. 

The mass and angular momentum loss rates for the 
different test cases are shown in the upper part of Ta- 
ble 2j For comparison and verification of these results, 
we have computed similar wind models for the solar case. 
The same numerical method described above was em- 
ployed for a solar magnetogram obtained during the last 
solar maximum (Carrington Rotation 1958) by the SoHO 
MDJfJ This solar simulation resulted in a mass loss rate 
of MfT, w 2 • 10- ■ 



-^^ Mq Yr- 



and an angular momentum 
loss rate of J « lO'^" g cm? s~^ with the use of base den- 
sity of 2 X 10^ cm~^ (same as Case A). These are similar 
to canonical solar values, as expected. 

Instead, the loss rates from AB Dor are significantly 
higher than solar. The angular momentum loss rate can 
be 2 orders of magnitude higher, simply due to the much 
more rapid rotation of AB Dor. For Case A, the mass loss 
rate for AB Dor is about a factor of 10 larger than the 
equivalent solar case. This is perhaps slightly surprising 
since all parameters other than the rotation rate and, to 
some extent, the surface field map are fairly similar to 
those of the active Sun chosen for the comparison. The 
most conspicuous difference is the factor of 50 in rotation 
rate and it is worth examining the influence of rotation 
alone in more detail. 

We repeated the AB Dor computational runs for Cases 
A-C for a rotation period of 25 d instead of 0.5 d, with 
all other aspects of the simulations remaining the same. 
The solutions for these runs are shown in Figure |4j The 
aziniuthal tangling of the coronal field that characterizes 
the 0.5 d period results in Figure |3] is, unsurprisingly, 
completely absent in this set of solutions and all field 
lines are essentially radial. In addition, more field lines 
are open in these solutions compared to the case with 
rapid rotation. Mass and angular momentum loss rates 
are listed at the bottom of Table [2] mass loss rates are 
typically a factor of ten lower than for the 0.5 d period 
results, and, for the Case A base density, are more similar 
to the solar value. 

4. DISCUSSION 



The MHD simulations of AB Dor presented here re- 
veals a coronal structure that is manifestly different from 
the well-studied solar corona. These differences are due 
to the different magnetic structure, which is mostly com- 
posed of high-latitude, large-scale regions of strong mag- 
netic field, and the rapid stellar rotation which induces 
azimuthal wrapping and tangling of magnetic field. This 
tangling cannot be obtained from the static, non-MHD, 
potential field extrapolation, which is generally useful 
only for studying the small closed loops near the surface 
where global effects are less important. We note in pass- 
ing that for stars with large-scale regions of strong mag- 
netic field, closed loops are probably significantly larger 
than in the solar case and the choice for the location of 
source surface location should be at greater radial dis- 
tance than the common use of Rss — 2.5 — 3.5i?*. 

The simulation results show that the mass loss and 
angular momentum loss rates increase with increasing 
coronal base density. The explanation for the former is 
trivial: introducing a greater mass source at the base will 
necessarily increase the mass flux through a closed sur- 
face around this source. The latter effect is more subtle 
and is due to the fact that J oc pu. When increasing the 
base density, the density drop with height decreases and 
the volume between the stellar surface and the Alfven 
surface is filled with more mass and as a result, more 
torque is being applied on the rotating star, thus increas- 
ing the angular momentum loss. 

For a given distribution of pu, the angular momentum 
loss rate is directly proportional VL^,, so it is not surpris- 
ing that J varies with the rotation rate. We find the 
mass loss rate also depends on the stellar rotation rate. 
This effect is not apparent if the rotation rate is not very 
high. However, the case of the extremely short rotation 
period of 0.5 d has a significant effect. The reason for 
this behavior can be found in the azimuthal tangling of 
the coronal field. When the rotation is slow, the global 
topology of the coronal field is radial (as seen in Fig- 
urel4]). In this case, the coronal density profile essentially 
drops like r~^. When strong rotation is present, the az- 
imuthal component of both the magnetic field and the 
flow become important and the radial component of the 
velocity is reduced. The increase in density in the slow 
wind case is greater than the decrease in speed, and as a 
result, the total value of pu increases. Another contribu- 
tor to the angular momentum loss increase with rotation 
is the shape of the Alfven surface for each case. 

Figures |5][7] show the shape of the Alfven surface for 
each test case with different rotation period. When the 
rotational period is small, the shape is modified to ac- 
count for the azimuthal component, and the surface is 
enlarged. In the case of fast rotation, the shape of the 
Alfven surface is modified and signs of the azimuthal 
component of the coronal field can be seen. It also seems 
like the size of the Alfven surface with faster rotation 
is slightly bigger. The simulations presented here show 
that the mass loss rate of AB Dor, and presumably also 
of other rapidly-rotating young late-type stars, is sub- 
stantially higher than the solar mass loss rate, and that it 
co uld be as high as 10~^^ — 1 0~^^ Mq Yr~^, as suggested 



by Wood et al. (2004 20051. These values however, are 



http : //sun . Stanford. edu/ 



strongly dependent on the assumed average coronal base 
density. It seems likely, based on measurements based 



on X-ray spectra presumably originating from plasma in 
closed loops (Section [2) that this is generally higher than 
the solar case, probaoly by an order of magnitude. For 
such a case, the predicted mass loss rate for AB Dor is 
about 100 times the solar rate. 

The mass and angular momentum loss rates found from 
our MHD models here are intriguing for the wider prob- 
lem of stellar rotational evolution. Models such as we 
present here could, in principle, be employed to map out 
theoretical AML as a function of stellar activity and ro- 
tation rate. While the general picture of stellar spin- 
down with age a s a result of wind-driven AML, emerged 



decades ago (e.g. Schatzmann|1962 Weber & Davis 1967 
Mestel 1968[ Skumanich 1972[ ), the details of situation. 



has proved somewhat complicated and rotation rate data 
amassed in the intervening years for late-type stars ex- 
hibits a complex dispersion over stellar age and mass. 

Faster rotation during stellar youth engenders greater 
magnetic activity through rotationally-powered dynamo 
action and the correlation of magnetic activity indica- 
tors such as chromospheric emission lines and coronal 
X-ray luminosity with rotation is well-established. Ob- 
servations of Lya absorption by the interactions of stellar 
winds with the surrounding ISM - the stellar equivalent 
of the heliopause - also indicate that stellar wind mass- 
loss rates are large r for younger and more active stars. 



Wood et al.| ( |2002[ ) estimate a relation M ex i-200±o.52^ 
based on combining inferred mass-loss rates with X-ray 
activity and observed X-ray activity vs stellar age. Their 
relation suggests that at very fast rotation rates, mass- 
loss should approach 1000 times the solar value, though 
they caution against the reliability of this extrapolation. 
Here, we find that plausible coronal base densities lead 
to mass loss rates of 100 times that of the present day 
Sun. 

Observations of rotation rates for stars in open clus- 
ters indicate that very young stars with ages of up to 
100 Myr or so are not as rapidly spun -down as would be 
expected based on extrapolation of the Skumanich ( 1972 1 
spin-down relation, and theoretical spin-down modelhng 
efforts have invoked a magnetic "saturation" that limits 
the AML for very short rotation period s (eg Ch aboye r 



et al. 1995 Krishnamurthi et al. 



1997 Barnes 



While our wind model tor a single star, such as Ab Dor 
presented here, cannot in itself be used to validate such 
a saturation approach, the general methodology does in 
principle allow for a more thorough exploration of the 
relevant parameter space to provide an MHD wind pre- 
diction of angular momentum loss as a function of stel- 
lar age. The requisite input parameters here would be 
the global coronal base density, and the magnetic field 
strength. 

The spindown time of AB Dor can be estimated as 
( |Ud-Doula eral][2009| : 



kM^Vt^Rl 



0.1 • X-^ir^ [yr] 



(7) 



where we used k — Q.l and write rA = XR^,. Based 
on Equation [7] and the calculated mass loss rates, the 
spindown time of AB Dor can range from 10^ — 10^^ yr 
depending on the particular case and on the value of X 
(which ranges from 5-10). A value of 10^ yr leads to a 
rotation of P = 0.5 d ■ e^ = 71 d after 5 billion years. 



However, Cohen et al. ( 2009[ ) have shown that the angu- 
lar momentum loss rate can be 3-4 times higher when the 
stellar magnetic field is dominated by strong polar spots, 
as appear to characterize young, fast-rotating stars. We 
expect the angular momentum loss rate to decrease with 
time as AB Dor becomes an established main-sequence 
star with spots at lower latitudes. Therefore, after five 
billion years a rotation similar to that of the Sun might be 
expected. In Equation [71 t is independent oi^^,. There- 
fore, for the same parameters but with different rotation 
rates, we have P1/P2 — 25 (i/0.5 d = 50 and J2I J\ = 50 
(the inverse of txJt-i). The ratios of the angular momen- 
tum loss rates for Cases A-C with different rotation rates 
are 60,68, and 80, respectively. This is consistent with 
the expected idealized value. 

In principle, we can propose a relation between the 
wind in our model and the Alfven radius as follows. Since 
at the Alfven radius, Usw = va we have: 



,{l/fs)=VA{rA) = 



y^AnpirA) 



VAoRl 



(8) 



where we assume that B(r) = _Bo(^*/?')^ and p(r) — 
pi^{R^/rY, with Bq, pq, and vao being the magnetic field, 
density, and Alfven speed at the flux tube base, respec- 
tively. From Equation [8] we get: 



rA 
R. 



VAO 



AVIsY 



(9) 



which could provide, in principle the location of the 
Alfven radius based on a known magnetic field distri- 
bution and surface density. The relation however, is not 
trivial due to the non-linear relation between the param- 
eters that define it. In reality, the radial functions of the 
magnetic field and density used in Equation IS] are not 
necessarily valid within the Alfven surface. 

5. SUMMARY AND CONCLUSIONS 

We have carried out a 3-dimensional global numerical 
MHD simulation of the corona of AB Dor, driven by a 
Zeeman-Doppler Image magnetic surface map. We stud- 
ied three test cases with different base density and we 
also compared the solutions with fast and slow stellar 
rotations. We find that the coronal structure of AB Dor 
is dominated by the azimuthal tangling of the coronal 
magnetic field as a result of rapid rotation. Based on 
the MHD solution, we calculate a realistic Alfven sur- 
face, which enables us to estimate the mass and angular 
momentum loss rates. Our main finding is that the mass 
loss rate is dependent on the value of the average coro- 
nal base density, as well as the coronal field and stellar 
wind topology which are affected by rapid stellar rota- 
tion. The total mass loss rates ranges between 10-500 
times the solar mass loss rate, while the angular momen- 
tum loss rate ranges between 15-30000 times the solar 
angular momentum loss rate. We demonstrate that the 
global coronal solution depends to some extent on the 
detailed properties of the coronal plasma. 

In addition to the uncertainty in stellar parameters, 
the surface maps used here to drive the model are not 
well defined as well. First, we interpolated the field in 
the "missing" part of the stellar surface, second the reso- 
lution of the maps in the regions where data is available 
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is not very high, and third, the interpretation of these 
maps is somehow debatable. In particular, large-regions 
of the map appear with strong magnetic field. One can 
ask whether these are really large-scale strong field re- 
gions, or whether it is more localized active region that 
is being smeared by the low resolution. The two scenarios 
should lead to a significantly different MHD solutions. 

Further modeling effort like the one presented here 
should focus on constraining stellar parameters such as 
the stellar wind speed for non solar-like stars. In addi- 
tion, a more consistent model to drive the stellar wind 
can help to generalize the model to stellar systems. 



We thank for an unknown referee for his/hers useful 
comments. OC is supported by SHINE through NSF 
ATM-0823592 grant, and by NASA-LWSTRT Grant 
NNG05GM44G. JJD and VLK were funded by NASA 
contract NAS8-39073 to the Chandra X-ray Center. Sim- 
ulation results were obtained using the Space Weather 
Modeling Framework, developed by the Center for Space 
Environment Modeling, at the University of Michigan 
with funding support from NASA ESS, NASA ESTO- 
CT, NSF KDI, and DoD MURI. 



REFERENCES 



Altschuler, M. D., Levine, R. H., Stix, M., & Harvey, J. 1977, Sol. 

Phys., 51, 345 
Altschuler, M. D., & Newkirk, G. 1969, Sol. Phys., 9, 131 
Arge, C. N., & Pizzo, V. J. 2000, J. Gcophys. Res., 105, 10,465 
Barnes, S. A. 2003, ApJ, 586, 464 
Budding, E., Erdem, A., Innis, J. L., Olah, K., & Slee, O. B. 

2009, Astronomischc Nachrichten, 330, 358 
Chaboyer, B., Deinarque, P., & Pinsonneault, M. H. 1995, ApJ, 

441, 876 
Cohen, O., Drake, J. J., Kashyap, V. L., & Gombosi, T. I. 2009, 

ApJ, 699, 1501 
Cohen, O., Drake, J. J., Kashyap, V. L., Korhonen, H., Elstner, 

D., & Gombosi, T. I. 2010, ArXiv e-prints 
Cohen, O., Sokolov, I. V., Roussev, I. I., Arge, C. N., Manchester, 

W. B., Gombosi, T. I., Frazin, R. A., Park, H., Butala, M. D., 

Kamalabadi, F., & Velli, M. 2007, ApJ, 645, L163 
Cohen, O., Sokolov, I. V., Roussev, I. I., & Gombosi, T. I. 2008, 

Journal of Geophysical Research (Space Physics), 113, 3104 
Collier Cameron, A., & Donati, J. 2002, MNRAS, 329, L23 
Cutispoto, G. 1998, A&AS, 127, 207 

Donati, J., & Collier Cameron, A. 1997, MNRAS, 291, 1 
Donati, J., Collier Cameron, A., Hussain, G. A. J., & Scmel, M. 

1999, MNRAS, 302, 437 
Donati, J., Morin, J., Delfosse, X., Forveille, T., Fares, R., 

Moutou, C, & Jardine, M. 2009, in American Institute of 

Physics Conference Series, Vol. 1094, American Institute of 

Physics Conference Series, ed. E. Stempels, 130-139 
Donati, J., & Semel, M. 1990, Sol. Phys., 128, 227 
Garcia- Alvarez, D., Drake, J. J., Lin, L., Kashyap, V. L., & Ball, 

B. 2005, ApJ, 621, 1009 
Grothues, H., Oestreicher, M. O., Gochermann, J., Tappert, C, 

Zaum, A., Brugger, H. R., & Schmidt-Kaler, T. 1997, A&AS, 

121, 247 
Giidel, M., Audard, M., Briggs, K., Haberl, F., Magee, H., 

Maggio, A., Mewe, R., Pallavicini, R., & Pye, J. 2001, A&A, 

365, L336 
Guirado, J. C, Reynolds, J. E., Lestrade, J., Preston, R. A., 

Jauncey, D. L., Jones, D. L., Tzioumis, A. K., Ferris, R. H., 

King, E. A., Lovell, J. E. J., McCuUoch, P. M., Johnston, K. J., 

Kingham, K. A., Martin, J. O., White, G. L., Jones, P. A., 

Arenou, F., Froeschle, M., Kovalevsky, J., Martin, C., 

Lindegren, L., & Soderhjelm, S. 1997, ApJ, 490, 835 
Hussain, G. A. J., Brickhouse, N. S., Dupree, A. K., Jardine, 

M. M., van Ballegooijcn, A. A., Hoogcrwerf, R., Collier 

Cameron, A., Donati, J., & Favata, F. 2005, ApJ, 621, 999 
Hussain, G. A. J., Jardine, M., Donati, J., Brickhouse, N. S., 

Dunstone, N. J., Wood, K., Dupree, A. K., Collier Cameron, 

A., & Favata, F. 2007, MNRAS, 377, 1488 
Hussain, G. A. J., van Ballegooijcn, A. A., Jardine, M., & Collier 

Cameron, A. 2002, ApJ, 575, 1078 
Janson, M., Brandner, W., Lenzen, R., Close, L., Nielsen, E., 

Hartung, M., Henning, T., & Bouy, H. 2007, A&A, 462, 615 
Jardine, M., Barnes, J. R., Donati, J., & Collier Cameron, A. 

1999, MNRAS, 305, L35 
Jardine, M., Collier Cameron, A., & Donati, J. 2002, MNRAS, 

333, 339 
Jiirvinen, S. P., Berdyugina, S. V., Tuominen, I., Cutispoto, G., 

& Bos, M. 2005, A&A, 432, 657 
Jeffers, S. V., Donati, J., & Collier Cameron, A. 2007, MNRAS, 

375, 567 



Krishnamurthi, A., Pinsonneault, M. H., Barnes, S., & Sofia, S. 

1997, ApJ, 480, 303 
Kuerster, M., Schmitt, J. H. M. M., & Cutispoto, G. 1994, A&A, 

289, 899 
Lim, J., White, S. M., Nelson, G. J., & Benz, A. O. 1994, ApJ, 

430, 332 
Lopez-Santiago, J., Monies, D., Crespo-Chacon, I., & 

Fernandez-Figueroa, M. J. 2006, ApJ, 643, 1160 

Luhman, K. L., Stauffer, J. R., & Mamajek, E. E. 2005, ApJ, 

628, L69 
Maggio, A., Pallavicini, R., Reale, F., & Tagliaferri, G. 2000, 

A&A, 356, 627 
Matranga, M., Mathioudakis, M., Kay, H. R. M., & Keenan, 

F. P. 2005, ApJ, 621, L125 

Mclvor, T., Jardine, M., Cameron, A. C, Wood, K., & Donati, J. 

2003, MNRAS, 345, 601 
Mestel, L. 1968, MNRAS, 138, 359 
Nielsen, E. L., Close, L. M., Guirado, J. C, Biller, B. A., Lenzen, 

R., Brandner, W., Hartung, M., & Lidman, C. 2005, 

Astronomischc Nachrichten, 326, 1033 
Pakull, M. W. 1981, A&A, 104, 33 
Phillips, J. L., Bame, S. J., Barnes, A., Barraclough, B. L., 

Feldman, W. C, Goldstein, B. E., Gosling, J. T., Hoogcveen, 

G. W., McComas, D. J., Neugebauer, M., & Suess, S. T. 1995, 
Gcophys. Res. Lett., 22, 3,301 

Pointer, G. R., Jardine, M., Collier Cameron, A., & Donati, J. 

2002, MNRAS, 330, 160 
Powell, K. G., Roe, P. L., Linde, T. J., Gombosi, T. I., & de 

Zeeuw, D. L. 1999, Journal of Computational Physics, 154, 284 
Sanz-Forcada, J., Maggio, A., & Micela, G. 2003, A&A, 408, 1087 
Schatzmann, E. 1962, Ann dAp, 25, 18 
Schmitt, J. H. M. M., Krautter, J., Appenzeller, I., Mandel, H., 

Wichmann, R., Barnstedt, J., Goelz, M., Grewing, M., Gringel, 

W., Haas, C, Hopfensitz, W., Kappelmann, N., & Kraemer, G. 

1997, A&A, 325, 249 
Skumanich, A. 1972, ApJ, 171, 565 
Toth, G., Sokolov, I. V., Gombosi, T. I., Chesney, D. R., Clauer, 

C. R., De Zeeuw, D. L., Hansen, K. C, Kane, K. J., 

Manchester, W. B., Ochmke, R. C, Powell, K. G., Ridley, 

A. J., Roussev, I. I., Stout, Q. F., Volberg, O., Wolf, R. A., 

Sazykin, S., Chan, A., Yu, B., & Kota, J. 2005, J. Geophys. 

Res., 110, 12,226 
Totten, T. L., Freeman, J. W., & Arya, S. 1995, J. Geophys. Res., 

100, 13 
— . 1996, J. Geophys. Res., 101, 15629 
Ud-Doula, A., Owocki, S. P., & Townsend, R. H. D. 2009, 

MNRAS, 392, 1022 
Vilhu, O., Muhh, P., Huovelin, J., Hakala, P., Rucinski, S. M., & 

Collier Cameron, A. 1998, AJ, 115, 1610 
Wang, Y.-M., & Sheeley, N. R. 1990, Astrophys. J., 355, 726 
Weber, E. J., & Davis, L. J. 1967, ApJ, 148, 217 
Wood, B. E., Miiller, H., Zank, G. P., Izmodenov, V. V., & 

Linsky, J. L. 2004, Advances in Space Research, 34, 66 
Wood, B. E., Miiller, H., Zank, G. P., & Linsky, J. L. 2002, ApJ, 

574, 412 
Wood, B. E., Miiller, H., Zank, G. P., Linsky, J. L., & Rcdfield, 

S. 2005, ApJ, 628, L143 
Zuckerman, B., Song, I., & BesseU, M. S. 2004, ApJ, 613, L65 



"0 


2- 10*, 10'-*, lOi'^ cm.-^ 


To 


5 MK 


-R* 


0.86/?Q 


M^, 


O.76Af0 


Prot 


0.5 d 



TABLE 1 
Adopted Properties of AB Dor 



TABLE 2 
Mass and Angular Momentum Loss Rates for AB Dor. 
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Fig. 1. — Longitude-latitude map of the original, incomplete ZDI based map of AB Dor (top-left), and the interpolated map used in the 
simulation displayed as a longitude-latitude map (top-right) and on two corresponding longitudinal spheres (bottom row). 




Fig. 2. — The global structure of the potential field extrapolation for AB Dor (left). The small sphere in the middle represents the stellar 
surface colored with grayscale contours of the radial field, while yellow lines represent the magnetic field lines. The outer white spherical 
shade represents the source surface located at r = 10-R*. Right panel shows a zoom close to the stellar surface. 




Fig. 3. — Numerical MHD solutions of the coronal magnetic field (top), wind speed (middle), and plasma density (bottom) for AB Dor for 
cases A-C (left to right with increasing base density). Top: global structure of the MHD solution with P* = 0.5 d. The sphere represents 
the stellar surface colored with contours of the radial magnetic field, the yellow lines represent the three-dimensional magnetic field lines, 
and the y = and 2 = planes are shown as transparent white shades. Middle: A side view of the top panel without the field lines. 
Bottom: A side view of the top panel with color contours of number density and zoomed around the star. The green shade represents an 
iso surface of n = 10* cm~^. Cases A-C are shown from left to right. 
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Fig. 4. — Same display as FigurelSlbut for the cases with P* = 25 d. 
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Fig. 5. — The Alfven surface for Case A from different view angle (left-right) colored with contours of the local value of M with P* 
(top) and with P^, = 25 d (bottom). White shades of the y = and z = plains are also shown. 
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Fig. 6. — Same display as Figurelslbut for Case B. 
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Fig. 7. — Same display as Figurelslbut for Case C. 



